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METHOD FOR CALCULATING THE ROLLING AND YAWING MOMENTS DUE TO 
ROLLING FOR UNSWEPT WINGS WITH OR WITHOUT FLAPS OR AILERONS 
BY USE OF NONLINEAR SECTION LIFT DATA^ 

By Aubbet P. Mabtina 


SUMMARY 


SYMBOLS 


The methods of NAG A Reports 866 and, 1090 have been 
applied to the calculation of the rolling-moment and yawing- 
moment coefficients due to rolling for unswept wings with or 
without flaps or ailerons. The methods are based on lifting- 
line theory and allow the use of nonlinear section lift data. 
The method presented herein permits calculations to be made 
somewhat beyond maximum lift for wings having no twist or 
continuous twist and employing airfoil sections which do not 
display large discontinuities in the lift curves. Oalculations 
can be made up to maximum lift for wings with discontinuous 
twist such as that produced by partial-span flaps or ailerons, 
or both. Two calculated examples are presented in simplified 
computing forms in order to illustrate the procedures involved. 

INTRODUCTION 

The calculation of the rolling-moment and yawing-moment 
coefficients due to rolling has received extensive treatment in 
the linear lift range; several of the better-known sources 
are references 1 to 3. Methods for making the calculations 
in the nonlinear lift range, however, are comparatively 
few (for example, refs. 4 to 6) and are based on the overall 
lift and drag of a wing. 

Methods for calculating wing characteristics in the non- 
linear lift range by using lifting-line theoiy and nonhnear 
section data (refs. 7 and 8) have given estimates which 
agree closely with results of wind-tunnel tests of nonrolling 
■wings, os seen in references 8 and 9. It was beheved, 
therefore, that these methods could be utilized in the cal- 
culation of the rolling derivatives and that such calculations 
might bo more accurate than those made by existing methods. 

Although the apphcation of the lifting-line method to a 
mUing wing in the nonlinear range is imphcitly contained in 
reference 7 and partially illustrated in reference 8, it is the 
purpose of the present report to outline the procedure of 
calculation by means of several illustrative examples. In 
addition, some new considerations regarding the apphcation 
of these methods to the calculation of the rolling character- 
istics are presented and discussed. Because of the supple- 
mentary nature of this report, the reader should be reasonably 
famihar with references 7 and 8. 


The term “section” as used herein denotes the section 
characteristics in three-dimensional flow. 
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■wing drag coefficient 
■wing rolling-moment coefficient 

coefficient of ■wing damping in roU, 

■wing lift coefficient 
■wing ya^wing-moment coefficient 
coefficient of ■wing yawing moment due to rolling, 

effective edge-velocity factor for symmetrical part 
of lift distribution, 


On 

Ct 

Or. 


E 


E' 


F 

K,K, 

R 

V 

Go 

b 

c 

C, 

c 

Cl 

CnSi 

P 


effective edge-velocity factor for antisymmetrical 

part of lift distribution, 1 

factor ■used in altering two-dimensional lift curves 
coefficients used in obtaining succeeding approxi- 
mations of lift distribution 
Beynolds ntunber 
velocity 

two-dimensional linear lift-curve slope 

■wing span 

local chord of wing 

root chord 

mean geometric chord, bjA 
section profile-drag coefficient 
two-dimensional profile-drag coefficient 
section lift coefficient 

additional section lift coefficient for (7i=1.0 
two-dimensional lift coefficient 
maximum section lift coefficient 
maximum two-dimensional lift coefficient 
section lift coefficient for a particular flap deflection 
angular robing velocity, radians/sec 


• Supersedes NAOA TN 2937, “Method for Oaloulatlng the RoUlng and Yawing Moments Due to Rolling for Unswept Wings With or Wlthont Flaps or Ailerons by Use of Nonlinear 
Section Lift Data" by Albert P. Martina, 19i3. 
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nmnber of intervals between points of calculation 
along span 

maxuniun thickness of wing section 
spanwise coordinate 
angle of attack, deg 

correction for induced angle of attack, deg 
effective angle of attack, deg 
induced angle of attack, deg 
angle of attack for zero lift, deg 
angle of attack for two-dimensional lift curves, deg 
angle of attack of root section, deg 
uncorrected induced angle of attack, deg 
induced angle-of-attack multiplier for asymmet- 
rical distribution 

magnitude of discontinuity in absolute and induced 
angles of attack, deg 

geometric angle of twist, negative if washout, deg 
approximate average angle of twist, 
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deg 


£p angle of twist due to rolling motion, deg 

6, angle of twist at wing tip, deg 

X taper ratio (ratio of tip chord to root chord) 

ij ratio of actual two-dimensional lift-curve slope to 

theoretical value of 

iin area multiplier for asymmetrical distributions 

v„ interpolation multiplier 

<r„ moment multiplier for asymmetrical distributions 


Superscript; 

* denotes value at end of flap or aileron 


CALCULATION PROCEDURE 

Inasmuch as complete theoretical developments of the 
method used herein are given in references 7 and 8, only the 
items pertinent to the purpose of the present report will be 
given. The method is based on lifting-line theory and uses 
the effective edge-velocity factors E and E' to correct for 
the effects of aspect ratio, following the concepts of reference 
10. Since E and E' as used herein are in themselves only 
approximations for elliptic wings, it should be realized that 
the results given by the present method at zero lift will not 
agree exactly with the values which would be obtained by 
using lifting-surface solutions for a given two-dimensional 
lift-curve slope. However, such differences are expected to 
be insignificant. 

The calculations for wings having no twist or continuous 
twist are discussed in the following section and are illustrated 
by TnAn.ng of example calculations for an rmtwisted wing 
designated as wing A. The calculations for wings with 
discontinuous twist require a slightly different procedme 
and therefore are discussed in an additional section. The 
latter calculations are illustrated by means of example cal- 
culations for a wing with partial-span flaps designated as 
wingB. 


WINGS WITH NO TWIST OH CONTINUODS TWIST 

The rolling motion of a wing produces an antisymmetric 
linear aerodynamic twist distribution which gives rise to 
the rolling derivatives. In the method of calculation used 
herein, the aerodynamic twist produced by the rolling motion 
is treated as a wing twist for purposes of obtaining the 
distribution of a<. Components of the inclined force vectors 
are then found, from which the rolling derivatives are 
determined. 

Use of section data. — ^The two-dimensional data at the 
Reynolds number appropriate to each wing station are simply 
plotted against an effective angle, a,=aoE- Either the 
section lift coefficients ci or the load coefficients cfi/b can 
be used, althou^ the use of the latter is believed to result 
in considerable savings in computing time. The load co- 
efficients were used in the example calculations for wing A 
and are presented in figure 1. -A single drag curve was used 
inasmuch as no account was taken of the spanwise variation 
in Reynolds number due to taper; the curve is presented in 
figure 2. The curves in figures 1 and 2 are based on data 
taken from reference 11 for the NAOA 65—006 section at a 
Reynolds number of 3.0X10®; however, the drag data were 



a*, deg 

Figdbb 1. — Section loading curves used in the example calculations for 

wing A. 
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Figurb 2. — Section drag curve used in the example calculations for 

wing A. 
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extrapolated to Mgher angles of attack on the basis of results 
in reference 12 for the NACA 64A006 section. 

Determination of the lift distribution. — The eflPeotive angle 
of attack at any station y is given by 


where 


and 


a,=a—oci—Aa, 

d) 

— a,+e+ep—at—Aoe, 

(2) 

180 py 
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(3) 


(4) 

= 2E' 

(5) 


The evaluation of equation (1) by the method of reference 
7 is one of successive approximations, in which a distribution 
of Cjc/6 is initially assumed. From this distribution and with 
the values of the multipliers dmt ^ven either in table I for 
r=10 or in reference 7 for r=20, a distribution of is 
obtained by using equation (4); equation (5) can then be 
evaluated. "With aU the components on the right-hand side 
of equation (2) determined, it is possible to find a,. Values 
of Cj corresponding to the values of are found from appro- 
priate section data, after which a check distribution of efijh 
can be obtained. If these check values of cfilh do not agree 
with those initially assumed, new values are assumed and the 
process is repeated imtil agreement is obtained. In the 


following examples, methods are indicated by which the 
differences between the check and initially assumed values 
are utilized to obtain the succeeding assumed values so that 
agreement is reached in a minimum number of approxima- 
tions. Equation (6) corrects the effective angle to accoimt 
for the 15' factor that is used with the antisymmetric com- 
ponent of loading when the section data are plotted against 

ag—aJE. 

Determination of the rolling-moment and yawing-moment 
coefficients. — ^After the distributions otcfifb and Ci^cjh are de- 
termined (distributions of Ca^c[h are obtained by using the 
plots of section data), components of these distributions 
along the x- and s-axes in the wind-axes system give the 
foUowing expressions: 

^=-^cos(«p — at)— ^sin(6p— a<) (6) 


C,C C|C , ■ . , . . \ 

SID (ep~afj 

For the angles usually encountered, 
cos (tp — «{)*« 1 


( 7 ) 


and 


sin («, — a<)«j^(«p— “<) 


so that equations (6) and (7) become 

CxC CiO V 


b 180 
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( 8 ) 


TABLE I.— INDUCED-ANGLE-OF-ATTACK MULTIPLIERS 
FOR ASYMMETRICAL LIFT DISTRIBUTIONS 
WITH r=10 1 
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It is easily seen that, for ep=0°, equations (8) and (9) reduce 
to the usual expressions for the section contributions to the 
lift and drag of a nonrolling wing inasmuch as the contribu- 
tion of Cog to the lift in equation (9) is usually negligible. 

Spanwise integrations of equations (8) and (9) yield the 
yawing-moment and rolling-moment coefficients as follows: 


‘^'180^^” b^ b 

[(■x)„~lfo (x)„ 


( 10 ) 

(10a) 

( 11 ) 

(11a) 


Values of the numerical integrating“multiphers cr„ are given 
in table II for r=10 and r=20 in addition to the area mul- 
tipliers r]„. The ijb, multiphers may be used to find Cl and 

Cd by similar relations, that is, Ol=A i-r-) Vm, aud so 

mol \ 0 


forth. 
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Application of the method. — ^In order to illustrate the use 
of the method, an example is presented for iving A at an 
angle of attack of 12°, rolling at such a rate that the ■wing-tip 
heKx angle p6/2F is 0.01 radian. The calculations are made 
for r=10. The pertinent dimensional data are given in 
table in and the lift distribution is calculated in table IV. 

The initially assumed load distribution in the linear lift 
range can he most rapidly obtained by using methods •with 
ivhich the reader is probably already familiar. For example, 
nearly exact initial assumptions can he made very quickly 
by proper use of the numerous lifting-surface solutions •which 


TABLE n.— WING-COEPjBTCIENT MTJLTIPLIEES FOR WINGS 
WITHOUT DISCONTINUOUS TWIST 
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are available in chart form and thus effect appreciable savings 
in computing time. For the sake of presenting a numerical 
procedure, ho^wever, the method of reference 13, ■which has 
been shown to give very good results in most cases, was used 
in this report. By combining the equation of reference 7 for 
the, initial approximation •with equation (33) of reference 8 
and equations (13) and (18) (modified for E') of reference 13, 
the foUo^wing expression is obtained for the initially assumed 
load distribution: 
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TABLE m.— GEOMETRIC CHARACTERISTICS OF WING A 

Taper ratio, X 0.60 Qeometrlc twist, «i, dec None 

Aspect ratio, A .4.00 Edgo-volodty tootor, E 1.118 
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Tip section NAOA 66-006 Wing Eeynolds number, I! S.OXIO* 
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TABLE rV.— CALCULATION OF LIFT DISTRIBUTION FOR WING A 
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where Cj(a) is the section lift coefficient for the geometric 
angle of attack in question. This loading is entered in col- 
umn 0 of table IV in normal order and in colmnn 0 in 
reverse order. The loa ding is entered in this manner in 
order to shorten the size of the computing form. The 
mechanics of computing are aU self-explanatory and it will 
be noted that the correction to the antisymmetric part of the 
loading is made in column @ as obtained by using equation 
(6). For the sake of brevity only the final calculations have 
been shoAvn in the table. The check load coefficients in 
column read from the section plots, will usually not agree 
with the assumed values for the first approximation. The 
process is repeated until agreement is obtained between 
the assumed and check values of cfijh. The manner used 
in determining the succeeding assumptions is dependent on 
a number of factors such as the linearity and the slopes of the 
section loading curves and on whether r=10 or r=20. 

Various methods for obtaining succeeding assumptions are 
presented in the appendix. When calculations are being 
made for more than one angle of attack, the first assumption 
for the second angle of attack can be based on the solution 
for the previous angle of attack by finding the value of Cjc/6 
corresponding to an angle 

“* 3 = “<1+ (“*j- “*i) ^ (13) 

in which it can be assumed that in the linear range. 

Once the values of or, have been determined for two angles of 
attack, plots of against a, can be made for each section. 

Values of cjc/6 corresponding to the extrapolated values of a, 

TABLE V.— CALCULATION OP THE ROLLING DERIVATIVES FOR WTNG A 
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wiU usually give a fairly accurate first assumption and thus 
minimize the amount of computing required. If only 
limited calculations are being made, it is recommended that 
any calculations in the nonlinear range be based on the 
results of a calculation in the linear range in a manner 
similar to that just described. In general, this procedure wiU 
eliminate a rather arduous solution since the load distribu- 
tions may change very rapidly in the nonlinear range. 

After the induced angles of attack and the lift distribution 
are determined, then the profile-drag distribution can be 
determined. Inasmuch as each section is assumed to be 
acting two-dimensionaUy, the section drag coefficients are 
obtained at the section lift coefficients or efilective angles a, 
for the proper values of Reynolds and Mach numbers. The 
calculations are carried out in columns 0 to 0 of table V. 
The calculations leading to the rolling derivatives are carried 
out in columns 0 to and at the bottom of the table. 

WINGS WITH DISCONTINUODS TWIST 

The discussion in this section is limited to the case where 
r=20; a similar method could be devised for r=10. It is 
emphasized that the method of altering the two-dimensional 
data, which is subsequently described, applies only up to and 
induding maximum lift and therefore precludes the calcula- 
tion of Oi^ and (7„p beyond maximum lift. 

Alteration of two-dimensional data. — The two-dimensional 
data to be used for wings with discontinuous twist must be 
altered in order to avoid a discontmuity in the spanwise 
distribution of maximum lift coefficient at the end of the 
flap or aileron since no such discontinuity exists in the 
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physical flow. The TnnxiTmiTn lift-coefficient values are 
altered (see ref. 8) by the relation 

where the spanwise variation of is given for several cases 
in figure 3. An illustration for deriving the F factors for 
any wing is given in figure 4. The quantity is the 

increment in cj^ at the discontinuity due to the deflection 
of the flap or aileron for the proper local Reynolds number. 
The values of Cj and oo are then altered according to the 
equations 


^t„„T 

(16) 


(16) 


The data altered in this manner are shown plotted in JSgure 
5 for wing B with 60-percent-span spht flaps deflected 60®. 
For purposes of comparison the rmaltered two-dimensional 
section data cross-plotted from reference 11 are also shown 
in figure 5. 

.The two-dimensional drag values as such are not altered, 
but the values of (Cdg)^ corresponding to the values of Ci^ 
are replotted against either Ci or o«. The drag data used in 
the example for wing B are shown in figure 6 plotted against 



a: 


Fiqube 3. — Factor for altering two-dimensional data for several wings 
having discontinuous twist. 


Cl. The data for the unflapped sections were taken from 
reference 11, cross-plotted to the values of B shown in 
figure 5. Since no drag data were available for the NAOA 
64r-210 section with spht flaps deflected, the data for the 
NAOA 23012 airfoil section (ref. 14) were used inasmuch as 
the lift curves were similar to those for the NAOA 64-210 
section up to maximum lift. Since the data of reference 14 
were for J?=3.5X10“, no accoimt was taken of the Reynolds 
number variation across the flapped portion of the span. 
The manner of alteration just described is necessarily arbi- 
trary and further experimental work may indicate a dif- 
ferent procedure; however, it should be recognized that the 
drag contributions depend on the differences in a, between 
the right and left wings and therefore are not critically 
dependent on the absolute values of the drag polars. 



Fiqueb 4. — Sohematio illustration of the calculation of the factor F 
used in altering the two-dimensional data for a wing having dis- 
continuous twist with the discontinuities located at ±2y*[b. 
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FiamiB 6. — Unaltered two-dimensional and altered section lift curves used in the example calculations for wing B. Two-dimensional curves cross 

plotted from data of reference 11. 
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Figubb 6. — Section drag curves used in the example calculations for wing B. 
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Determination of the lift distrihntion. — ^For wings with, dis- 
continuous twist the induced angle of attack computed by 
means of equation (4) is modified by a term which is pro- 
portional to the magnitude of the discontinuity and acts as a 
correction factor to account for the inability of a limited 
trigonometric series to represent adequately the lift dis- 
tribution of a vtong with partial-span flaps or ailerons de- 
flected. This induced angle, according to reference 8, is 

“<*=««»+ 

where is the uncorrected induced angle of attack given 
by equation (4), 5 is the magnitude of the discontinuity, and 
values of the correction factors Oej5 taken from reference 8 
axe given in table VI. These correction factors depend only 
on the spanwise position of the discontinuity and apply to a 
deflected flap or aileron extending from the point of dis- 
continuity to the right wing tip. Values for any spanwise 
extent of the deflection of the flap or aileron, or both, can be 
obtained by combinations of the values presented. For 
example, if the flap deflection is symmetrical and extends over 
60 percent of the wing span, the values of a^/5 would be 


obtained by subtracting the values of aJS corresponding to 

2v* 2v* 

-|-=0.6 from those which correspond to -|-=— 0.6. The 

tabular values of aJS for ^or the flap side of the 

discontinuity therefore the corresponding section 

lift curves should be used with these values. The value of 5 
is obtained as shown in flgiue 5. 

If the discontinuity is located at a station othei’ than 

^—cos^i the v^ues of (cic/b)* are interpolated from the 

calculated values of Cjc/6 as follows: 


values of the interpolation multiphers are given in table 
YU for various spanwise positions of the end of the flap. 
Multiphers for positions of 2y*/b other than those tabulated 
can be calculated from the relations given in reference 8. 

Application of the method. — The method is apphed to 
wing B with 60-percent-span spht flaps deflected 60° at 
«,=10° and rolling at such a rate that the wing-tip helix 
angle pbj2V is 0.01 radian. The geometric characteristics 


TABLE VT.— ANGLE CORRECTION ^ AT STATION ^ DEE TO DISCONTINEITY IN INDUCED ANGLE OF ATTACK AT 

3 0 

STATION ^ 

0 


\ 

\6 

2 (f\ 

-0.9877 

-a95U 

-a 8910 

-0.8090 

-a 7071 

-a 5878 

-0.4640 

-a3090 

-a 1664 

0 

aii« 

0.3000 

0.4M0 

0.5878 

a 7071 

0.8090 

0.8910 

a 9611 

0.0877 

-a 0877 

0.4340 

-a 0262 

a0029 

-00013 

O.OOH 

-a 0003 

a 0001 

-a 0001 

0 

QHH 

0 



0 


M 


0 


-.0611 

.0768 

.4666 

0324 

.0039 

-.0018 

.0006 

-.0004 

.0002 

-.0002 





0 




0 


-.9000 

-.0213 

.1U4 

.2648 

-.0036 

.0022 

.0006 

.0002 

.0003 

.0001 

H 

0 

.0001 

^BBB 





0 


-.8910 

-.0134 

.0612 

.4786 

-.0358 

.0044 

-.0021 

.0007 

-.0005 

.0002 

B txiy 

■iMTiT 

-.0001 

B 

B^^3 



KpH 

0 


-.8090 

.0078 

-.OICB 

.0661 

.4849 

-.0379 

.0048 

-.0023 

.0008 

-.0006 

B 

BIm!1/ 

.0001 

B txo 

B^^3 

B 

.0001 


0 

-.0001 

-.8000 

-.0035 

-.0053 

.0072 

-.3362 

-.0657 

.0068 

-.0051 

.0012 

—.0014 

B 

Hfrm 

.0002 

B LjJA 

B^^S 

B ^2 

.0001 

Bl» »;vl 


-.0003 

-.7071 

-.0028 

.0058 

-.0087 

.0517 

.4890 

-.0396 

.0051 

-.0026 

.0008 

B 1x0 

■tTfi 

-.0003 

B Wij 

B^^3 

B LjW 

-. 0001 

.0001 


.0001 

-.7000 

0016 

-.0011 

-.0060 

.0164 

-.3802 

-.0639 

.0066 

-.0051 

.0012 

B uff 

■fixi! 

-.0007 

B LjW 

B^^S 

B ujj 

-.0003 

.0001 

E)^2 

.0001 

-.6000 

.0063 

-.0036 

.0120 

-.0U7 

.0876 

.3154 

-.0071 

.0030 

.0009 

B 

HtiWi 

.0001 

B IJJA 

BBB 

B ^2 


.0002 

^QBi 

.0003 

-.3878 

.0021 

-.0021 

.0049 

-.0079 

.0496 

.4920 

-.4408 

.0054 

-.0027 

B ixff 

B 

.0003 

B ^2 

B 

iB iSSL 

.0001 

-.0001 

B| ^ 

-.0001 

-.6000 

-.0038 

.0061 

-.0069 

.0124 

-.0183 

.0630 

.0964 

.0133 

—.0025 

B 

B 1 IlK 

.0016 

B jSjj 

B ^2 

B ^2 

.0006 

-.0003 

m 

-.0003 

-.4640 

-.0010 

.0016 

0017 

.0043 

-.0074 

.0480 

.4944 

-.0419 

.0057 

B 

B nin 

-.0008 

B ^2 

B ^2 

B ffiy 

-.0002 

.0001 

^9 iViy 

.0001 

-.«00 

.0023 

-.0039 

.0031 

-.0064 

.0056 

-.0190 

-.0666 

-.0413 

.0146 

B 

B ikI 

-.0026 

B 

B L[J£ 

B u2 

-.0009 

.0007 

^9 rm 

.0000 

-.3090 

.0003 

-.0007 

.0013 

-. 0016 

.0039 

-.0070 

.0467 

.4964 

-.0428 

B 

B i.i'ii 

.0011 

B SS 

B lW 

B Wji 

.0002 

-.0002 

iViv 

-.0003 

-.3000 

-.0007 

-.0003 

-.0008 

-.0008 

-.0006 

-.0049 

.0162 

-.3786 

—.0690 

B rw 

B ms 

.0016 

B hW 

B ^2 

B loi 

.0004 

-.0007 

^B m 

-.0000 

-.2000 

-. 0017 

.0029 

-.0022 

.0046 

-. 0042 

.0107 

-. 0163 

.0H3 

.1839 

.0128 

B IIIK 

.OOU 

B 

B lW 

B ^2 

.0013 

-.0005 

^B fit 

-.0005 

-. 1664 

-.0004 

.0006 

-.0006 

.0011 

-. 0013 

.0036 

-.0067 

.0467 

.4983 

-.0438 

Brr!i 

HfTigTp 

B Mjf 

B 

B KMi 

-.0004 

.0003 

^9 m 

.0003 

-.1000 

.0012 

-.0020 

.0015 

-.0027 

.0023 

-.0052 

.0013 

-.0170 

—.0776 

—.0491 

.0167 

-.0088 

B 

B j2i 

B Wf 

.0020 

.0013 

^9 m 

.0010 

0 

.0004 

-.0003 

.0005 

-.0005 

.0010 

-.0012 

.0034 

-.0064 

.0447 


-. 0447 

.0064 

B 

B 

B 

.0006 

-.0006 

^B m 

-.0004 

.1000 

-. 0010 

.0016 

.0013 

-.0020 

-.0019 

.0034 

-. 0040 

.0088 

-.0167 

.0401 

.0776 

.0170 

B 

B LuV 

B lWj 

.0027 

-.0016 

^B jjW 

-.0013 

.1664 

-.0002 

.0003 

-.0003 

.0004 

—.0004 

.0009 

-. 0011 

.0032 

—.0061 

.0438 

mLMM 

-.0467 

B ^2 

B ^2 

B 

-.0011 

.0006 

^9 iVm 

.0004 

.2000 

.0006 

-.0010 

.0006 

-. 0013 

.0007 

-.0020 

.0012 

-.0041 

.0018 

-.0128 

-.1339 

-.0713 

B ilTs 

B 

B 

-. 0046 


^9 nwj 

.0017 

.3000 

.0006 

-.0003 

.0007 

-.0004 

.0011 

-.0006 

.0021 

-. 0015 

.0063 


.0600 

.3786 

-. 0162 

B uMj 

B ^2 


.0008 

^B iVr< 

.0007 

.3090 

.0002 

-.0002 

.0002 

-.0002 

.0004 

-.0004 

.0008 

— .OOU 

.0031 

BfiW? 

.0428 


-.0467 

B rirai 

B 22 

.0016 

-.0013 

^B fW 

-.0008 

.4000 

-.0006 

.0008 

-.0007 

.0009 

-.0009 

.0014 

-. 0016 

.0026 

-.0034 

KS 

-. 0146 

.0413 

.06M 

.0190 

B jjw 

.0064 

-.0031 

^B mu 


.4540 

-.0001 

.0002 

-. 0001 

.0002 

-.0003 

.0003 

-.0004 

.0003 

-.0010 

BEjS 

B [52 

.0410 


-.0480 

.0074 

-.0043 

.0017 

^9 mfi 

.0010 

.6000 

.0003 

-.0005 

.0003 

-.0006 

.0004 

-.0009 

.0006 

-. 0016 

.0012 

B 22! 

B ijWj 

-.0133 



.0183 

-.0124 

.0059 

^9 nM 

.0038 

.6878 

.0001 

-.0001 

.0001 

-.0001 

.0002 

-.0002 

.0003 

-.0003 

,0007 

B ^2! 

B LjS 

-.0054 

.0406 

B^^l 

-.0496 

TiTiil 

-. 0W9 

^B LjW 

-.0021 

.6000 

-.0002 

0 

-.0002 

0 

-.0002 

0 

-.0003 

-. 0001 

-.0006 

B ixS 

B 


.0071 

-.3164 

-.0876 

.0117 

-. 0120 

^B iTtVi 


.7000 

-. 0001 

.0002 

-.0001 

.0003 

-.0001 

.0004 

-.0002 

.0007 

-.0004 

B 

B 

.0061 

-.0066 



-.0164 

.0060 

^B fnt 

.0016 

.7071 

-.0001 

.0001 

-.0001 

.0001 

-.0001 

.0001 

-.0001 

.0003 

-.0003 

B rfo 

B cj5j 

.0026 

-.0051 

.0396 

.6110 

-. 0517 

.0087 

^9 iVjs 

.0028 

.8000 

.0002 

-.0001 

.0002 

-.0001 

.0002 

-. 0001 

.0003 

-.0002 

.0006 

B ^2 

B LJjfc 

-. 0012 

.0051 


.0667 

.3362 

-. 0072 

^B ixi 

.0036 

.8090 

.0001 

0 

.0001 

-. 0001 

.0001 

-.0001 

.0001 

-.0001 

.0002 

B ixiy 

B 

-.0008 



.0379 

.6151 

-.0651 

^B 

-.0078 

.8910 

0 

0 

0 

0 

0 

.0001 

-.0001 

.0001 

-.0001 

B iTi? 

B ^2 

.0005 


.0021 


.0368 

.6215 

-. 0612 

.0134 

.9000 

0 

0 

0 

0 

0 

-. 0001 

0 

-.0001 

0 

B ojT 

B 

-.0003 



B^23 


-.2648 

-. 1114 

.0213 

.0611 

0 

0 

0 

0 

0 

0 

0 

0 

.0001 

B ixii 

B ixi? 

-.0002 

be^si 


BE2B 





.9877 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Hfl 

BH 

.0001 

gmi 

Biii 

bm 

.0013 
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of the wing are given in table 'VJJUL and the calculation of the 
spanwise lift distribution is given in table IX. The deter- 
mination of the initial assumption for the lift distribution is 
not given bereia inasmuch as the procedure is fully illustrated 
for this wing in reference 8 for the symmetrical case. For 
the asymmetrical case the approximate antisymmetric 
loading given by equation (18) of reference 13 (modified for 
E') 


CiC 

b ^2 _ 


--H.273. 




AE'+i 


should be added to the assumed symmetrical loading given 
by the procedure used in the example of reference 8. Any 
suitable method for obtaining the assumed loading may be 
used, however. 


TABLE Vn.— INTERPOLATION MULTIPLIERS FOR 

OBTAINING VALUES OF AT THE END OF THE FLAP 


2tf\ 6 

0.2000 

0.3000 

0.4000 

0.6000 

0.6000 

0.7000 

OAOOO 

0.9000 

0 

-a 2970 

-0.1839 

0 

0 

0 

0 

0 

0 

.1664 

LOGOS 

.3203 

-.1940 

0 

0 

0 

0 

0 

.3090 

.4209 

.8967 

.6203 

-.1274 

0 

0 

0 

0 

.4640 

1307 

-.0321 

.6946 

.7663 

0477 

-.0337 

0- 

0 

.6878 

0 

0 

-.1208 

.4786 

.9112 

.1234 

-.0307 

0 

.7071 

0 

0 

0 

-.U76 

.1878 

.9133 

.1634 

0 

.8090 

0 

0 

0 

0 

06U 

-.0330 

.9277 

-.0431 

.8910 

0 

0 

0 

0 

0 

0 

-.0444 

.8901 

.0611 

0 

0 

0 

0 

0 

0 

0 

.1951 

.0877 

0 

0 

0 

0 

0 

0 

0 

-. 0131 


The values of in column @ of table IX are those for 
0 — 


-=—0.6 from table VI: 


whereas the values of t-t in 

5+ 


colunm @ are the negative of those for -y-=0.6 and were 

obtained by the procedure outlined in the previous section. 
The values of 5 at ± 2y*Jb axe found from the flapped and 
unflapped section curves (fig. 5) at values of ci* correspond- 
ing to the values of (cjc/6) *. After these approximate values 
for S are determined, the- check span load distribution can 
then be calculated. The check v^ues of (cic{b)* are then 
interpolated from equation (18). When convergence is 
established (i.e., column @ minus column © ia zero), the 
values of at ±2y*/& can be found. (The method given for 
case I of the appendix was used to obtain convergence be- 
tween the assumed and check values of cic/b, columns © and 
respectively.) 

The drag distribution is calculated in columns © to © of 
table X. Note that two values are given at -^=±0.60. 

2i/ 

The values at -y= TO. 60 TO coirespond to the flapped sides. 

The rolling-moment and yawing-moment components are 
calculated in columns © to ® and the coefficients and 
derivatives are calculated at the bottom of table X. The 
multiph^ (T„ used in table X are tabulated in table XT. 
If the values of Cr, and Qd are desired, the n-umerical inte- 
grating area multipliers of reference 8 may be used. 


TABLE Vm.— GEOMETRIC CHARACTERISTICS OP WING B 


Taper ratio, X 0. 400 

Aspect ratio, A 9. 021 

Span, 6, a 16.000 

Root chord, e., K 2. 381 

Wing Reynolds nomber, B 4. 44X10“ 


Root section NAOA 61-210 

Tip section NAOA 64-210 

Grometrlo twist, ii, dM —2. 00 
Edge-velocity toctor, B 1-024 

Edge-velocity factor, Bf L 094 



c I 
^ 1 

R 

a 1687 

6.98X108 

.1438 , 

6.42 

.1293 

4.87 

.1166 

4.35 

.1027 

3.87 

.0914 

a44 

.0817 

ao8 

.0739 

2.78 

.0681 

2.67 

.0590 

2.22 

.1016 

1 

as3 



For tapered wings with stralght-Une elements from root to construction tip: 
e, b ti c/c, 

(Alter -values of eje, near tip to allow for (Use value of etc. before rounding tip.) 
rounding.) 


0 

0 

00006 

-14 

0 

.0016 

-.30 

0 

.0031 

-.60 

-.0477 

-.0063 

-.73 

.9112 

.3164 

-.98 

.1876 

.0878 

-L28 

-.05U 

-.0117 

-1.63 

0 

.0122 

-1.77 

0 

-.0030 

-L94 

0 

.0066 
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TABLE EX.— CALCULATION OF LIFT DISTRIBUTION FOR WING B 


Fo-e-span split flaps deflected 60°; <t,=10.00°; ^-0.01 radian! 



2»i I /CK 



® 

® 



MoltipUers 

16 

15 

14 

13 

12 

-.809 

-.707 

-.588 

-.454 

-.309 



1 

® 

cte 


T 

Reverse 

am 

(•) 

of® 


«ai637 


.. (t): 


-.951 0 


-.89l| 0 


-.0511 .0336 -26.374 


-166.985 

0 

-7.019 

463.533 

-122.749 

0 




-97.524 202.571 -71.139 0 -6.391 


-12.604] 0 I -81.392 177.054 - 64.735 


-71.29& 160.761] -60.725 


-64.817] 150.611] -58.514 


-60.768 145.025 


-58.533 


-6.530 


-2.19 


-6.192 .214 


.2431 


.891 


,951 


.988 0 




-1.068 


-.604 0 I -.981 


-.297 0 -.452 



- (?): 
G); 


a+= 12.07 



.951 

.891 

V ’ 

2 

3 


.707 

.588 

5 

6 


.309 

.156 

8 

9 
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& 

Q 

® 

@ 

® 

o 

® 

B 


a . 

a , 

®X3- 

® XS + 

(or+-*-f- 

a< 

Aa« 

G) 

«• 

@ 

Cl 

(Sectloii 

data) 

1 

1 

1 

o 

8.91 

—4.72 

-a 02 

13.66 


B 

0.1637 


7.29 


1.64 



aos 

7.40 

4.60 


3.01 

.467 


.0276 



-.m 

7.68 

Z21 



.723 

.0631 

.0402 



.16 


wm 



.896 


.0661 

-.0117 

0 

-.14 

8.28 

-.29 

-.02 

8.69 

1.026 


.0837 

.0378 


l.OS 

8.61 

-1.34 



Mwm 


.1053 

mm 


■EE 

8.&^ 

K! 

mm. 



10^ 

.1627 

-.0071 


-.08 

o.:m 

4.33 

-.01 

4.60 

1.826 

.1156 

.2100 


.0001 

.04 

9.6^ 

4.27 

-.01 


L869 

.1293 

.2404 

,0009 

.0008 


0.77 

4.78 

0 

4.09 


.1438 

.2636 

.0003 


.01 


6.89 

0 

4.U 

1.743 

.1687 

.2766 

.0006 

.0009 


9.96 

4.83 

0 


1.813 

WPS 

.2660 

.0001 


.04 

9.87 

4.34 

.01 


1.880 


.2431 

.0003 

-.0071 

-.03 

9.78 

4.44 

.01 

6.31 

1.868 

1^ 

.2146 

lElillill 

.31U 

3.81 

0.61 

t68 


191 

L628 

WIffI 

.1670 



108 

■BE 


m 

BS 

1.207 

1^ 

.1103 

■Qm 


-. 14 

9.31 

j|B|Q 


9M 

1.087 

1^ 

.0888 

,0002 


.18 

&08 

.98 

m 

7.98 

.961 



0 

-.0036 

-.04 

8.77 

2.68 


fljQR 


Rj^l 

.0630 

HQ 

.0063 

.03 

&63 


■ 

3.46 

.514 


.0303 

¥-0 

■ 

7.43 

■ 

■ 

■ 

.1016 

.1682 ' 

f+o 


■ 
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TABLE X— CALCULATION OP THE ROLLING DERIVATIVES FOR WING B 


[0.6-spaa split flaps deflected 60“; a,— 10.00°; ^-0.01 radton] 


® 

® 

® 

® 

® 

® 

® 

® 

® 

® 


® 

He 

b 

Cl 

(table IX) 

"o 

( flg . 6 ) 

e 

b 

b 

®x® 

0 . 673X ® 

at . 

(table EQ 

<•— at 

S-o 

b ^ 

CtC 

b 

(table IX) 

—X 

(table XD 

-a 983 

a467 

a0064 

a0690 

a00032 

- 0,67 

4.60 

- 6.07 

- 0.0016 

a0278 

-a 1410 

- 0.00809 

-.951 

.723 

.0092 

. 068f 

.00063 

-.64 

2.21 

-a 76 

-.0017 

.0492 

-.1368 

-.00709 

-.691 

.895 

.0103 

.0739 

.00076 

—.61 

.66 

- L 17 

-.0009 

.0662 

-.0774 

-.02118 

-.809 

1.025 

.0115 

.0317 

.00094 

-.46 

-.29 

-. 17 

-.0002 

.0836 

-.0143 

-.01272 

-.707 

1. 162 

.0123 

.0914 

■EH 

-.41 

- L34 

.94 

.0011 

.1053 

-.0990 

-. 0M98 

-. 600-0 

L613 

.0147 

.1016 

.00149 

—.34 

- 4.72 

4.38 

.0066 

.1637 

.6733 

-.00534 

-. 000+0 

1.513 

.1731 

.1016 

.01769 

-.34 

7.29 

- 7.65 

-. 1342 

.1637 

-41737 

-.00088 

-.638 

1.681 

.1730 

.1027 

.01777 

-.34 

6.54 

- 6.88 

-.1222 

.1627 

- L1194 

-.00716 

-.464 

L826 

.1730 

. U66 

.01908 

—.26 

4.33 

- 4.59 

-.0917 

.2109 

-.9680 

-.02118 

-.309 

1.869 

.1730 

.1293 

.02237 

-.18 

4.27 

- 4.46 

-.0995 

.2404 

- L0C93 

-.00709 

-.166 

1.833 

.1780 

.1433 

.02488 

-.09 

4.78 

-487 

1211 

.2635 

-4 2833 

-.00809 

0 

1.713 

.1731 

.1687 

.02747 

0 

5w80 

- 6.89 

-. 1618 

.2767 

- 1.6298 

0 

.166 


.1730 

.1438 

.02488 

.09 

4.83 

- 8.74 

-.0930 

.2650 

-.9911 

.00806 

.309 

LSSO 

.1730 

.1283 

.02237 

.18 

4.34 

- 4.16 

-.0931 

.2431 

-.40113 

.00700 

. 4M 

L858 - 

.1730 

.1166 

.01998 

.26 

4.44 

— 4.18 

-.0836 

.2146 

-.8970 

.02118 

.683 

L628 

.1730 

.1027 

.01777 

.34 

6.68 

- 6u34 

-. 1126 

.1671 

-40594 

.00710 

. 600-0 

L657 

.1730 

.1016 

.01758 

.34 

7.43 

- 7.09 

-.1246 

.1682 

-41216 

.00088 

. 600+0 

1.657 

.0155 

.1016 

.00168 

.84 

- 4.64 

4.98 

.0078 

.1583 

.7878 

.00534 

.707 

L207 

.0132 

.0914 

.00121 

.41 

-LM 

L55 

.0019 

.1103 

.ino 

.03498 

.309 

1.037 • 

.0128 

.0817 

.Ml® 

.46 

-.03 

.48 

.0005 

.0888 

.0420 

.01373 

.891 

.961 

-OUO 

.0739 

.00081 

.51 

.08 

-.47 

-.0004 

.0710 

-.0334 

.03118 

.951 

.778 

.0096 

.0681 

.00065 

.64 

aos 

- 2.14 

-.0014 

.0531 

-.1136 

.00709 

.933 

.614 

.0070 

.0690 

.ooou 

.57 

6.15 

—4 68 

-. 0019 

.0301 

1379 

.00809 

c.«-a2®x[®+^®] — o.ooin 

_ Cl - a00471 . 

2V 

C. ®3 — 0.00133 

' pb om 
2V 


DISCUSSION 

Lack of eitiher experiment -wing-alone rolling data or 
suitable two-dimensional section data prevents the makmg 
of exact correlations of. this method with experiment in the 
vicinity of maximum lift. It was possible, however, to 
compare the wing-alone calculations for -wing A -with experi- 
mental -wing-body results, siuce sufficient section data 
existed to allow the calculations to be carried slightly beyond 
maximum lift. The comparisons are presented in figure 7. 
Agreement is considered to be good when the differences be- 
tween the conditions of the calculations and tests are con- 
sidered. The failure of the present method to predict the 
increase in (7jp for 5° <a<8° may be inherent iu the method 
since similar experimental trends have been observed for 
other -wings of this plan form. The differences between the 
calculated and experimental variations of G^j, in the same 


angle-of-attack range could be partly associated -with the 
pre-viously noted iacrease m Oij, and partly due to body- 
interference effects as sho-wn in reference 16. 

Calculations which included the body-induced anglo-of- 
attack distribution on the -wing were made for ^ving A up to 
10° angle of attack. The indications were that this com- 
ponent of body interference had negligible effects on the 
calcffiated variations presented m figine 7. Additional 
calculations for -wing A using a section drag of C|tan a (no 
leading-edge suction) yielded a Onp variation that agreed 
quite closely -with the experimental variation in the low 
angle-of-attack range; such a (7,^ variation would indicate 
the possible existence of early separation in the experimental 
'results, a likely possibility considering the low test Reynolds 
number and the section thickness ratio. The differences in 
the high lift range, particularly -with regard to the (7„p vaiia- 
tions, probably can be aU ascribed to body-interference 
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TABLE XI— ROLLING- AND YAWING-MOMENT-COEFPICIENT MULTIPLIERS FOR WINGS WITH DISCONTINUOUS 

TWIST AND r=20 AND HAVING DISCONTINUITIES LOCATED AT 

0 

[Values far poslUve ^ are shown; for negative reverse signs of all multipliers] 


\ 

\ b 
^ \ 

0 

a 1604 

aso90 

a 4510 

a6878 

aTOTl 

aS090 

a 8910 

a 9511 

a9877 

(•) 


0.1000 


aooTOO 

a 00314 

a 02118 

a 01218 

a02618 

0.01246 

a 03118 

a 00769 

a00809 

a00033 

-a 00067 

.iGei 

0 


• OHM 

.01986 

.01246 

.02618 

.01246 

.02118 

.00769 

.00809 

.00202 

.00253 

.i!000 

0 

.00766 

.00898 

.02061 

.0X216 

.02618 

.01246 

.02US 

.00769 

.00809 

-.00255 

.00233 

.3000 

■■ 

.00735 

.00420 

.02U8 

.01216 

.02618 

.01245 

.02118 

.00769 

.00809 

.00340 

.00033 

.3090 

■■ 

.00800 


.02118 

.01246 

.02618 

.01246 

.02U8 

.00769 

.00809 

.00385 

.00385 

.4000 


.00794 

.00816 

.01831 

.01318 

.02618 

.01246 

.02118 

.00769 

.00809 

.00389 

-.00203 

.4640 

0 

.00768 

.01164 


.01867 

.0^64 

.01246 

.03118 

.00769 

.00809 

.00662 

.00662 

.6000 

0 

.00809 

.00812 

.01883 

.01383 

.02660 

.01245 

.02118 

.00769 

.00809 

-.00378 

.00491 

.6878 

0 

.00809 

.00769 

.02118 


. .02818 

.01246 

.02U8 

.00769 

.00809 

.00623 

.00623 

.6000 

0 

.00809 

.00760 

.02118 

.oono 

.02498 

,01372 

.02118 

.00769 

.00309 

.00088 

.00534 

.7000 

0 

.00809 

.00769 

.02006 

.01770 

.06188 

.01268 

.03U8 

.00769 

.00809 

.00771 

-.04770 

.7071 

0 

.00800 

.00769 

.01938 

.01807 


.01887 

.01986 

.00769 

.00809 

.00818 

.00818 

,8000 


.00809 

.00769 

.02118 

.01372 

.02498 

.00716 

.02118 

.00789 

.00809 

.00534 

.00083 

.8000 

■■ 

.00809 

.00769 

.03118 

.01246 

.02618 


.02118 

.00769 

.00809 

.00623 

.00623 

.S910 


.00809 

.00769 

.02118 

.01246 

.03454 

.01868 


.0U64 

.00759 

.00662 

.00662 

.0000 

mm 

.00809 

.00769 

.02U8 

.01246 

.02618 

.01280 

.02968 

.01034 

.00773 

-.01673 

.00567 

.0511 


.00809 

.00769 

.02118 

.01246 

.02618 

.01246 

.02118 


.00809 

.00385 

.00385 

.0877 

m 

.00809 

.00769 

,02U8 

.01246 

.02618 

.om5 

.01988 

.0U54 



.00253 

.00304 


•ForneBBtfro ~j this bMoma— — 0^ - — +0. 

** For nesattvo j tHli become*— —0. 


Method Conflg. Airfoil 


o Exp. ■ 

Cole. Ret 5 

Cole. Ret 6. 


Wing-body NACA 65A006 


Cole. Present Wing A. NACA 65-006 
report 


R 


0.7X10® 

14 

3.0 

13 



Figubb 7. — The calculated and experimental variations of and 
with angle of attack for two wings having aspect ratios of 4.0 and 
taper ratios of 0.60. Experimental data from the 6-foot-diameter 
rolling-flow test section of the Langley stability tunnel. 


effects since the differences axe similar in trend to those of 
reference 15. 

Also shown in figure 7 are the variations of Oj^, and Cnp 
calculated by means of references 5 and 6, respectively, in 
which use is made of the experimental wing-body lift and 
drag curves. It appears that both the present method and 
that of reference 5 give reasonable estimates of the variation 
of Otp. The present method, however, appears to give a 
more realistic picture of the variation of in the high 
angle-of-attack range than that of reference 6. The edge- 
suction effects which are appHed in reference 6 were negligible 
for this aspect ratio and therefore do not account for the 
differences shown. 

Several items of interest were observed during the course 
of the calculations for wing A. For example, the contribu- 
tions of the section lift and profile drag to the rolling deriva- 
tives could be separated as shown in figure 8. The profile- 
drag component of is seen to be opposite and nearly two 
to three times that due to the lift at high angles of attack. 
The profile-drag contribution to the damping in roU, on the 
other hand, is either zero or negligible even at the higher 
angles of attack. 

The method presented herein is primarily intended for use 
in the nonlinear range where some flow separation is present, 
and the subsequent discussion briefly covers some considera- 
tions which limit the use of this method. The condition 
that must be essentially fulfilled is in keeping with the basic 
assumption that all sections operate two-dimensionally or 
nearly so. Therefore any separation that is present must 
not give rise to large amounts of spanwise flow which can 
cause a complete redistribution of lift and thus invalidate 
the basic assumption. Lack of sufficient experimental data 
on load distributions beyond maximum lift prevents the 
formulation of positive limits concerning these phenomena, 
but, on the basis of airfoil characteristics and on observed 
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PiaxTRE 8. — The calculated rolling derivativea for wing A showing the 
lift and profile-drag components. K=3.0X 10*. 

stalling behavior, regions most likely to be amenable to 
calculation can at least be classified. 

It is believed that departures from the two-dimensional- 
flow assumption will not seriously affect the calculation of the 
rolling derivatives for wings which incorporate airfoil sections 
that exhibit gradual changes in lift-curve slope beyond 
maximum lift. On the other hand, calculations beyond 
maximum lift for wings which incorporate sections having 
large discontinuities in the lift curves are beheved to have 
little significance because with such airfoils there is no known 


way of predicting either the extent of the initial stall or the 
influence of the stalled area on the section characteristics of 
neighboring sections. In addition, such wings display a 
tendency toward asymmetrical stall imder no-roll conditions 
which leads to rolling and yawing moments. The deflection 
of trailing-edge flaps generally produces abrupt lift-ciuwe 
peaks at maximum lift on all but the thin airfoils. Such 
conditions coupled with the inability to treat the discon- 
tinuous-twist cases beyond maximum lift, as previously 
discussed, limit the feasibility of calculations beyond maxi- 
mum lift to a very limited number of cases with full-span 
flaps. 

Although calculations may appear to be feasible on the 
basis of airfoil-section stalling characteristics, there is an 
additional consideration which may limit the extent of the 
calculations. This limit, referred to in reference 16 as the 
stability limit, is 

that is, an increment in the angle of attack Aa cannot result 
in a decrease in the induced angle by an amount greater 
than the iucrement Aa. An idea of the meaning of this 
limit can be obtained from the calculations for wing A. The 
maximum negative value in this case was —0.76 compared 
with —1.0 given by equation (19). 

Although one has an apparent choice of performing calcu- 
lations with either r=10 or r=20 for the unflapped case, it 
has been fmmd from experience that calculations beyond 
maximum lift with r=20 rarely are required except perhaps 
at very high aspect ratios. If such calculations are neces- 
sary, however, it is suggested that the results from calcula- 
tions with r=10 be used as the initial approximations for 
the r=20 solution, a procedure which will generally shorten 
computing time. 

Langley Aeronautical Laboratory, 

National Abvisory Committee for Aeronautics, 
Langley Field, Va., January S9, 195S. 



APPENDIX 


METHODS FOR OBTAINING SUCCEEDING APPROXIMATIONS FOR THE SPANWISE LOAD DISTRIBUTION FROM AN 

INITIALLY ASSUMED LOAD DISTRIBUTION 


METHODS FOB USE WITH r-20 

These methods for obtaining succeeding approximations 
for the load distribution have been extensively used and 
produce convergence in reasonably few approximations. 

Case I: positive and linear.- — ^For a positive and linear 

lift-curve slope, the succeeding approximation, denoted by 
the superscript 1, is given by the equation 



0 2 4 6 8 10 12 


V 

Figtob 9. — Coefficients used to obtain succeeding approximations. 

r=20. 


where 




in which the increments A | 


are the differences between 


the check values and the assumed values (column @ minus 
column ®, table IX), and K and are constants for any 
particular wing. Values of K and Ki taken from reference 
8 are presented in figure 9 as functions of AJEji). Values of 
K{ for i greater than 3 are small enough to be considered 
negligible. 

The number of terms of equation (A2) needed for any 
particular approximation depends upon the convergence of 
the assumptions; fewer terms are needed when the differences 

A small or when positive differences nearly cancel 

negative differences. Equation (A2) applied to wing B at 

fUz=0 ^^=0.688^ becomes 

A' {3.2(®-®),-k((g)-(D)*-l-(@-©)7+ 

0.4 [(@-®)4+(@-®)8] + 0.3 [(@-®)3+(@-®)9] } 
For m— 1 and 2, equation (A2) is expanded as 

(f ) (f ) +(K-XSi (f) + 

aA(f)+X.A(f)J (A3) 

(f H (f)+^ (f) + 

since 

For m=18 and 19, the values of A' 

similar in form to the values of A' 
respectively. For wing B, equation (A3) would be 

A' p.2-0.4) (@-®)i-h 

(1.0-0.3)(@-®)2-h0.4(@-®),-l-0.3(@-®)J 
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whereas equation (A4) would become 

A' (x) t(l-0-0-3)(@-(D)i+3.2(@-©h+ 

1.0(@— ®)3+0.4(@— ©)4+0.3(@— ®)5] 

Use of the K factors has been fomid to establish con- 
vergence within three approximations when the mitial 
assumption is reasonably dose. .. . 

Case n : positive and nonlinear. — ^Although the factors 

in figure 9 were derived for a linear lift-curve slope that is 
constant across the span, they can be used in the nonlinear 
range. An estimation is made of* the wing lift-cuurve slope 
at the angle of attack in question, from which a value of -n is 
obtained. The K factors are then obtained at the proper 
value of AEJt]. In this way it is possible to use these factors 
for values of jj as low as 0.3 for A=3.0. From the trends 
of the variations, it is seen that considerable extrapolation 
is permissible at the higher values of AElrj. 

Case HE: negative and linear or nonlinear. — ^For a 

negative lift-curve slope, either linear or nonlinear, the 
procedure is as follows: 

(1) First obtain values of the span load distribution with 
fewer significant figures than desired. 

(2) Find A station, and work with these 

values directly. 

(3) Adjust the load distribution at those stations where 

the largest values of A (^x) adjusting the loadings 

at m and n? ± 1 stations to obtain approximate convergence 

at the wth station. If large values of A ^x) 

adjacent stations, then adjust the loading at only one of the 
stations, although with some practice the adjustments 
required at two such stations become quite obvious. 

(4) As the values of A made smaller, the nmnber 

of significant figures in the solution can be increased to the 
desired amount. 

(5) Keep all values of A ^x) positive or negative, 
if possible, for easier mental manipulation. 

METHODS FOB USE WITH r=10 

Although factors similar to the K’s given in case I could 
be derived for the corresponding condition with r=10, it is 


equally facile to use the method described under case IH 
for all cases with r= 10. 
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